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The  dynamics  of  solvation  is  of  great  importance  for  many  aspects  of 
solution  chemistry,  such  as  spectroscopy,  reaction  kinetics,  and  dynamic  light 
scattering.  It  is  not  surprising,  then,  that  a  number  of  excellent  discussions 
exist  on  the  subject^ .  The  availability  of  experimental  data  for  fast 
reactions  in  solution  has  prompted,  very  recently,  theoretical  research  or.  the 
effects  of  the  discrete  structure  of  the  solvent  on  the  relaxation  of  internal 
modes  in  the  solutes,  either  ionic  or  neutral  species.  In  particular,  ve  should 
mention  the  work  of  Wolynes^^  and  of  Rips,  Klafter  and  Jortner  who 

discuss  the  relaxation  of  the  shift  of  an  effective  hydration  sphere  diameter, 
using  the  mean  spherical  approximation  (MSA)  as  the  theory. 

In  the  present  discussion  we  want  to  study  a  different  approach  to  this 
problem;  namely,  the  electric  microfield  distribution  function  and  its  dynamic 
behavior. 

The  electric  microfield  distribution  was  first  studied  by  Holtsmark^"^ , 
who  studied  the  weak  coupling  regime.  We  are  interested  in  the  formulation  of 
the  microfield  distribution  function 


W(0  -  W(e.t)  t_Q  -  <  6((  -  Kq(0)> 


-  |  dJRV  dp 


■v  e'^H  6((  -  Eq) 


(2) 


where  we  are  using  canonical  ensemble  averages  of  the  3-dimensional  Dirac  delta. 

V  CF 

function  $(<■).  Our  system  has  v  particles  of  coordinates  *1*  *2  -  “n-R  and  m 

Y 

momenta  —  p^  -  p  .  The  inverse  temperature  /9  -  1/kgT  where  kg  is 
Boltzmann's  constant  and  T  is  the  absolute  temperature.  The  Hamiltonian  of  the  ~~ 
system  is  H.  We  single  out  a  particle  0,  located  at  the  origin  -  p^  -  0. 


The  electrostatic  field  produced  by  a  particle  j ,  of  coo 


Codes 

^tes  Rj  ,  Pj  is  Eqj  ,  id/ or 


□  □ 
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Furthermore 


vv 

j-i 


'oj  -  •  v  W 


where  <ciqj<Rj)  is  the  electrostatic  potential  at  the  origin  produced  by  particle 

j .  Normally  the  evaluation  of  the  electric  microfield  would  be  done  by 

expanding  Eqn.  (2),  which  entails  the  calculation  of  all  v-particle  distribution 

functions .  The  problem  could  also  be  reformulated  in  terms  of  a  path  integral , 

as  has  been  done  very  recently^.  But  a  more  practical  approach  was  proposed  by 
17  18 

Morita  ,  and  more  explicitly  by  Iglesias  .  In  this  formalism  the  Fourier 
transform  of  W(£) , 

A(A)  -  f  d(  eiX  *  •  V(0  (5) 

is  introduced.  It  can  be  shown  that  the  calculation  of  A(A)  can  be  performed 

using  the  pair  correlation  functions  of  a  modified  hamiltonian  H(A).  Explicit 

results  have  been  obtained  in  the  MSA  for  the  one  component  plasma  on  a 

19 

neutralizing  background 

The  static  microfield  problem  for  ionic  and  polar  solutions  was  discussed 
20 

recently  .  In  the  MSA,  and  for  hard  ions  and  spherical  solvent  molecules,  the 
general  result  is  that  the  microfield  distribution  at  the  center  of  a  spherical 
test  charge  is  Gaussian,  of  the  form 

3  1 — 3 —  -  2  i! 

U  «>  -  J  2^  •  2  "2  <6> 


where  m2  -  <(  >  is  the  mean  square  average  of  the  field  at  the  testing  particle. 
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In  several  cases  of  interest  the  value  of  ■ 2  is  known  explicitly,  such  as 
in  the  For  a  Mixture  of  ions  of  charge  z^e  (|e|  -  electron  charge), 

density  p ^  and  hard  core  diaaeter  o A  one  obtains 

®2  “  +  \  pl  zi  g0i(a0i)  (7) 


where  ggi^Oi^  *s  C^e  P*^r  c°ntact  distribution  function,  and  z^  is  the  charge 
of  the  test  particle  located  at  the  origin.  In  the  Onsager  limit,  for  which  the 
Debye  screening  length  is  zero,  the  following  simple  and  exact  result  holds: 


-  <U  .> 
/3  fi-1 


(8) 


where  <U^>  is  the  Onsager  self  energy  of  a  dipole  ft.  In  our  case  ft  -  1.  For 
the  primitive  model  electrolyte  (8)  reduces  to  (7),  since  the  MSA  second  moment 
is  exact  in  the  Onsager  limit. 

The  first  observation  is  th<  t  the  second  moment  m^  is  also  the  fluctuation 
of  the  electrostatic  energy,  <E  •  K>.  Then  the  natural  extension  of  the 
microfield  idea  to  non-equilibrium  situations  is  to  study  the  quantity 

L(t)  -  <E(0)  •  K(t)>  (9) 

where  E(t)  is  given  by  eqns.  (3,  4)  with  coordinates  Rj(t)  (t  "  time); 
therefore,  L(t)  is  the  fluctuation  (and,  by  the  Nyquist  theorm  also  the 
dissipation  function)  of  the  solvation  energy .  Introducing  the  quantity 

W(^  e2,  t)  -  <6^x  -  E0(0)]  -  E0(t)]>  (10) 

we  see  that,  after  a  short  calculation 


L(t) 


w(eL,  t2,  t)  ^  t2 


<E(0)  •  E(t)> 


(ID 
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we  imagine  now  a  situation  in  which  our  teat  particle  has  undergone  some  change, 
in  which  its  c,  *rga,  or  better,  its  charge  distribution  has  changed.  Assuming 
that  the  change  is  permanent  the  system  will  initially  be  in  a  state  in  which 
the  caicrofleld  distribution  is  Vq(£).  At  the  end  of  the  process  the  aicrofield 
distribution  will  be  tf  (£).  The  Implication  is  that 

40 


and 


tf<*l  <2*  t)|t<_0  "  V*l’  ♦>  *  «2> 

w<e1  ?2.  t)|  -  w0(^,  4)  v„  <e2) 


Furthermore,  ,  t)  must  also  obey  the  normalization 


(12) 


(13) 


^2  W(^l  ^2* 


w(f1) 


(14) 


at  any  time  (t).  These  conditions  imply  that  the  conditional  microfield 
probability  distribution  is  given  by 


V(«Q.  t,  t) 

K(£,  tjf0)  -  0) 


(15) 


and  K(f,  tjf^)  relaxes  from  a  i((  •  (q)  function  to  the  Gaussian  distribution 


K(f.  «U0) 


(16) 


The  time  dependence  of  K(£,t|£g)  is  complicated  because  of  the  complexity 
of  the  system.  It  is  therefore  reasonable  to  assume  that  it  is  given  by  some 
Gaussian  random  process  (this  does  not  mean  that  £  is  a  Gaussian  random 
variable)  which  is  characterized  by  a  propagator  and  a  source  term.  It  can  be 
verified  that  the  boundary  conditions  (12-16)  are  satisfied  if 
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K<€,  tUQ)  -  dtL  C^)  r<t  -  tx) 


G< tL)  - 


(4xDt^) ' 


4Dt. 

•  1  .  A* 


r(t  -  t.)  -  — — ■— —  cos 

1  /  8*D  8)r  «2D 


W  |  -JL_  .  I 

4  2Dm„t  4 


2C  2D°2tl 


3  (19) 


The  cine  dependence  enters  through  the  diffusion  constant  D,  which  characterizes 
the  random  relaxation  process. 

An  alternative  approach  is  to  utilize  the  Fourier  series  representation  of 

22 

a  Gaussian  randon  process  ((  is  the  zero-mean  Gaussian  random  variable)  to 
derive  a  convenient  and  provocative  form  for  the  Joint  probability  distribution 

p(e.  t|e0): 


P(f,  tUn)  -  K  <«,  tjc0)  exp 


-  3  c 


where  the  transition  probability  K,  or  equivalently,  the  conditional  probability 
density  K  fixed),  is  given  by 


K(£.  tUn)  - 


2m»2(l  -  M  (t) ) 


exp  - 


€2  +  H2(t)  •  2  e0-  (  K(t) 

1  -  M2(t) 


Here  M(t)  is  the  normalized  microfield  autocorrelation  function 
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M(t) 


<((0 )  •  ((0)> 


(22) 


where  <  >  denotes  an  equilibrium  ensemble  average.  Note  that  Eqn.  (21)  is  a 

dynamic  contact  transformation  in  which  the  operator  K  acts  on  the  Maxwellian 

distribution  to  produce  P((,  1 1  (q) ,  Berne,  Pechukas,  and  Harp  have  shown  that 

23 

this  result  can  also  be  obtained  from  an  information  theoretic  approach 
Therefore,  if  (  is  a  Gaussian  random  variable, 


tU0) 


lV2 


■  U  - 


M2(t)) 


exp 


3 

e2  ♦  *2  -  2  (  .  (0  M(o 

2“2 

1  -  M2(t) 

t,  4 

• 

Furthermore,  if  M(t)  decays  as  a  simple  exponential 
M(t)  -  exp  (at)  (a  >  0) 


(22) 


(23) 


Chen  we  recover  an  Omstein-Uhlenbeck  (OU)  process,  which  is  "essentially" 

(allow  for  linear  transformations  of  (  and  t)  the  only  process  which  is 

24 

stationary,  Gaussian,  and  Markovian  .  It  follows  that,  if  in  a  microfield 
model  such  as  the  ion-dipole  MSA,  one  generates  an  M(t)  which  decays  as  a  sum  of 
simple  exponentials,  or  if  in  more  elaborate  models  one  encounters  sums  of 
oscillatory  exponentials,  then,  given  that  the  (-fluctuations  are 
time-stationary,  they  cannot  be  both  Gaussian  and  Markovian,  and  in  these  cases 
it  is  the  former  condition  which  does  not  hold. 

More  exotic  forms  forms  for  M(t)  such,  as  the  algebraic  decay 


M(t)  -  t  large,  (a  >  1)  (24) 

t 

or  stretched  exponential 

M(t)  -  exp(-t^), 


t  large,  (0  <  0  <  1) 


(25) 
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are  not  anticipated  for  electric  field  fluctuations  in  staple,  non-glassy  media. 
Algebraic  decay  lntplles  that  the  Harkov  assumption  has  broken  down,  while 
stretched  exponential  behavior  indicates  the  (more  serious)  apparent  violation 
of  time  stationarlty. 

As  a  staple  exaaple  of  non-Gaussian  (  fluctuations,  one  might  consider  a 
so-called  "dichotomic"  Harkov  process  (DHP)  in  which  the  aicrofleld  (  takes  on 
the  values  ±  (q  with  equal  probability .  In  this  case 

K«.  C|f0)  *  J  [l  ('«)]  i(  t(o 

+  |  |l  -  exp  (-ot)j  6 £  ^  (26) 

where  a  is  twice  the  transition  frequency  between  states  and  6  is  the  Kronecker 
delta  function.  Just  as  for  an  OU  process,  the  autocorrelation  function  of  a 
DHP  is  given  by 

H(t)  -  exp  (-at)  (27) 

Therefore,  if  only  two  point  correlation  functions  are  considered,  very  little 
can  be  said  about  the  stochastic  process  which  generates  aicrofleld 
fluctuations,  and  it  is  for  this  reason  that  we  have  emphasized  the  importance 
of  P(f,  t|£Q)  and  the  conditional  probability  density  K(^,  t|£0). 
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